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Electrically driven single-photon emitting devices have immediate applications in quantum 
cryptography, quantum computation and single-photon metrology. Mature device fabrication 
protocols and the recent observations of single defect systems with quantum functionalities 
make silicon carbide (SiC) an ideal material to build such devices. Here, we demonstrate 
the fabrication of bright single photon emitting diodes. The electrically driven emitters 
display fully polarized output, superior photon statistics (with a count rate of >300 kHz), 
and stability in both continuous and pulsed modes, all at room temperature. The atomic 
origin of the single photon source is proposed. These results provide a foundation for the 
large scale integration of single photon sources into a broad range of applications, such as 
quantum cryptography or linear optics quantum computing. 


With breakthroughs in SiC growth technolo¬ 
gies and its excellent thermal, mechanical and 
physical properties^, SiC has become an out¬ 
standing wide band gap semiconductor for both 
emerging optical and electronic applications. 
Well-developed CMOS compatible processing 
protocols exist for SiC providing the availability 
of broadband photonic optical cavities^ and, 
with the possibility of controlled doping of 
both p- and n-type impurities®, commercialized 
electronic devices. The compound nature and 
extensive polytypism of SiC results in a broad 
array of defects some of which have promising 
applications in quantum computing, cryptog¬ 
raphy and metrology. Optically driven single 
photon emission in 4H-SiC has been demon¬ 
strated for the carbon-vacancy carbon-antisite 
CyCsi pair^, the silicon vacancy (Vsi)^’® and 
the divacancy (VcVg;)', and in 3C-SiC for the 
CyCsi pair®. Other single photon emitters in 
wide band-gap materials such as ZnO® may also 
display similar characteristics. 

To significantly enhance the practicality 
of single photon emission, it is desirable to 


integrate single photon emitters into a device 
to allow electrical pumping. Electrically driven 
quantum dots have been demonstrated®®^®^ 
but these mostly operate at low temperatures. 
The possibility of electrical excitation at room 
temperature of the nitrogen vacancy defect in 
diamond, NV®, has also been shown®®’®6 
present there are severe limitations on the types 
of devices that can be engineered with this 
material, namely, the inability to create a high 
active n-type concentration by ion implanta¬ 
tion. In this work, single photon emitters in the 
visible spectral range are integrated into 4H-SiC 
p+n and n+p junction diodes. Results with 
6H-SiC appear in the supplementary material. 
The fabrication of these single photon emitting 
diodes is achieved by applying well developed 
CMOS compatible processes including pho¬ 
tolithography, ion implantation and annealing. 
Photon anti-bunching measurements indicate 
the quantum nature of the single defect. The 
emitters have high emission rates, are fully 
polarized and operate with high stability at 
room temperature, characteristics which are 
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FIG. 1. Single photon sources in bulk 4H-SiC and fully fabricated p+n junction diodes, a, 

Schematic of the confocal set-up used to characterize the single photon emitters. It includes a Hanbury, 
Brown-Twiss interferometer with two single photon avalanche detectors (D1 and D2) connected to a corre¬ 
lation card (CC). The dichroic mirror (DM) was removed when in used in EL mode. A partial schematic 
of the device consisting of a p’*' top contact is also shown. Three floating guard rings encircle the central 
contact to decrease the electric held at the main contacted junction, b, IV-curve of the diode. In this 
device, features due to shunt resistance, tunnelling current, diffusion current and series resistance are in¬ 
dicated by the letters a-d, respectively, c, EL map of the edge region of a device, d. Room temperature 
EL spectrum showing the source of the background (Di line) and the single photon emitter with a ZPL at 
745 nm. e, Background corrected anti-bunching traces with = 0) < 0.1 indicating excellent single 

photon emission characteristics. 


beneficial for quantum cryptography protocols 
and linear optics quantum computing^^. 

Results 

Single photon source and device integra¬ 
tion. SiC p'’'n junction diodes were formed 
in 4H SiC n-type epi-layers by implanting alu¬ 
minum while the substrate was held at 800° C. To 
activate the dopants a post-implantation anneal 
of 1600-1800°C was employed. During this an¬ 
neal a carbon capping layer was used to prevent 
Si sublimation and step bunching formation^®. 
No irradiation step was performed to create sin¬ 
gle defects after annealing. Fig. la shows the 
device schematic (bottom) together with the 
confocal microscope used to observe the mate¬ 
rial’s luminescence. The dopant profile is de¬ 
rived from SRIM calculations which are pre¬ 
sented in the Supplementary Information. A 


typical current-volt age (IV) curve from our fab¬ 
ricated diodes is displayed in Fig. lb and shows 
the expected rectifying behavior of a pn junction 
diode with a reverse leakage current of around 
1 X 10“^*^ A/cm“^. In forward bias, the de¬ 
vice displays clear features that are dominated 
by the series resistance which limits the injected 
current. Electroluminescence (EL) is also ob¬ 
served when driven above the threshold voltage 
of V ~ 3.1 V. Eurther fabrication and measure¬ 
ment set-up details are provided in the Methods 
section below. 

The EL arises when minority carriers are in¬ 
jected into the n-type region and radiatively re¬ 
combine as previously observed via deep-level 
defect states^®. Eig. Ic shows an EL map of the 
device under an applied forward bias of 10 V. 
The central bright spot is a single defect with 
a count rate of 250 kHz and does not display 
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FIG. 2. Single defect PL in 4H-SiC at 80 K. a, Room temperature confocal map excited with 532 nm. 
b,c, Single defect spectra taken at 80 K under excitation with 532 nm and 633 nm for two different defects, 
respectively. TO and LO indicate Raman lines. The ZPL widths are 0.5 nm for green and 0.8 nm for red 
excitation, d, Histogram created from more than 50 spectra indicating the number of lines per 5 nm interval. 


any blinking behaviour. The defects are located 
near the SiC epilayer surface between the top 
aluminum contact and the floating guard rings 
and not in the implanted regions. 

Figure Id shows a typical EL spectrum taken 
at a single defect site. The spectrum displays 
two main components. Firstly, a signal arising 
from the Di center with a weak zero phonon 
line (ZPL) at 427 nm and a phonon-side band 
(PSB) with a peak intensity around 520 nm is 
observed. This often appears as background lu¬ 
minescence in the vicinity of the implanted re¬ 
gion and is mostly formed as a result of the im¬ 
plantation induced damage, surviving the high 
temperature activation anneals. The Di line is 
associated with a Si antisite (Sic)^°. Secondly, 
in this example, ZPL emission at a wavelength 
of 750 nm in Fig. Id from the diffraction limited 
bright spot appears. To construct the EL map of 
single photon emitters shown in Fig. Ic, emission 
was collected through a 560 nm long-pass filter 
to suppress the Di line. Fig. le shows an exam¬ 
ple of the background corrected^^ anti-bunching 
trace with = 0) < 0.1 confirming excellent 

single photon emission characteristics. 

While the Di line is relatively well known, a 
wide variation of ZPL positions and line shapes 
was observed for single defects (see Supplemen¬ 
tary Information). To aid in the identifica¬ 
tion of these defects we studied their formation 
in a broad range of SiC materials with differ¬ 


ent doping densities and processing conditions. 
For these studies photoluminescence (PL) mea¬ 
surements on high-purity semi-insulating (HPSI) 
SiC from CREE was performed and enabled a 
greater number of stable defects to be studied 
at once. Additionally, when laser excitation was 
used with the single photon emitting diodes, the 
single defects could not be clearly characterized 
as their photo-stability decreased and a signifi¬ 
cant background luminescence appeared. 

The as-purchased HPSI sample showed only 
unstable emission as observed previously'^. As 
mentioned above, the defects of interest in this 
work form and are stabilized via the high tem¬ 
perature anneal only. As with the EL emission 
described in Fig 1, over 90% of these centers are 
completely photo stable. Further, it was found 
that if the surface roughness is allowed to de¬ 
grade during a high temperature anneal without 
the carbon capping layerthe photo-stability 
and defect density greatly decreased. The de¬ 
fects are thus very sensitive to the state of the 
surface and detailed confocal maps confirm that 
they exist in a narrow band close to the surface 
(see Supplementary Information). 

A typical PL confocal map of the HPSI SiC is 
shown in Fig. 2a. The defect density was found 
to vary slightly for different samples typically in 
the 2 — 8 X 10^ cm“^ range. Interestingly, elec¬ 
tron irradiation of 2 MeV electrons to a fluence 
of 1 X 10^^ cm“^ did not appear to significantly 
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increase the defect density or the type of defects 
observed after the high temperature annealing. 
By lowering the temperature to 80 K, most de¬ 
fect spectra develop sharp ZPL, with line widths 
ranging from 0.3-2 nm and could be character¬ 
ized further. Fig. 2b and c show two typical 
single defect spectra taken at 80 K for excita¬ 
tion wavelengths of 532 nm and 633 nm, respec¬ 
tively. Room temperature spectra are presented 
in the Supplementary Information. Phonon as¬ 
sisted transitions are still present as can be seen 
by the elevated background signal. The trans¬ 
verse optic (TO) and longitudinal optic (LO) 
phonon Raman modes are indicated. We per¬ 
formed spectral measurements of more than 50 
defects and created a histogram of the recorded 
ZPL peak positions (Fig. 2d). The histogram 
clearly shows the variety of ZPL over the whole 
red spectral region. We did not observe any 
groups of lines as might be expected for differ¬ 
ent charge states, defects on inequivalent lat¬ 
tice sites or multiple defect types. Instead the 
lines are continuously distributed over a 200 nm 
range. The large spread is peculiar and will be 
addressed in the discussion below. In the follow¬ 
ing sections we examine several other properties 
shared by all observed single photon emitters in 
these materials to give further scope to the dis¬ 
cussion. 

Single photon emitter dynamics. Figure 3a 
and 3b show the background corrected anti¬ 
bunching traces of two different single defects in 
EL (p''“n device) and PL (HPSI sample) mode 
that have similar spectral characteristics, respec¬ 
tively. A typical signal-to-background ratio is > 
0.85. Photon bunching around r = 10 ns is evi¬ 
dent in EL and especially in PL modes which is 
related to a dark metastable state. EL and PL 
modes differ most strikingly around r = 0 where 
the EL anti-bunching exhibits a plateau char¬ 
acteristic of a double pump process^®’^^. Such 
a process may involve, for example, two step 
carrier capturing process. The equation that 
is used to fit the PL data is a double expo¬ 
nential which is commonly used for three level 
systems^^. A standard analysis based on such 
a system yields a relatively fast excited state 
lifetime of = 3.3 ± 0.3 ns with transitions 
to a dark state every 10 cycles. All other cen- 



FIG. 3. Electrically driven single photon emit¬ 
ter dynamics. Background corrected a, EL and 
b, PL traces for two spectrally similar single pho¬ 
ton emitters in 4H-SiC (solid line fit described in the 
text), c, EL traces for different currents offset 

vertically for clarity. 


ters studied show very similar transition rates 
= 2 — 5 ns). 

To account for the plateau that is evident in 
the EL data, a third exponential term is required 
and yields an excellent fit. Such an empirical 
relation has a form that may be reflective of a 
four level electronic system^®’^^. Fig. 3c illus¬ 
trates the (r = 0) evolution under increasing 
injected current. Due to the faster pump rate, 
the = 0) sharpens and there is a subtle 

increase in the photon bunching. All time con¬ 
stants associated with the fit increase linearly 
with injected current. 

The intensity of the single defect emission in 
EL varies strongly from defect to defect with 
many being well over 300 kHz. One contributing 
factor to the intensity variation is the distance 
from the p^n or n'’'p junction which determines 
the current density at the defect location. A 
typical EL intensity-current curve is shown in 
Eig. 4a. After subtracting the linear background 
that originates from the Di center emission (dot¬ 
ted line), we obtain a linear relationship. Note, 
that the initial non-linearity is due to the un¬ 
known current density at the defect site and does 
not indicate a saturating behavior. The detec- 
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FIG. 4. Intensity-excitation power curves for single photon emitters in 4H-SiC. Intensity curve 
for a single defect as a function of a, current (EL) and b, laser power (PL) corrected for the back-ground 
(dotted line). The raw signal is indicated by the dashed line. 


tor count rate reached 350 kHz at 14 mA without 
signs of saturation. This intensity is more than 
twice as bright as the saturation count rate of an 
optically driven single nitrogen vacancy (NV“) 
in diamond investigated using the same experi¬ 
mental set-up, and more than 6 times brighter 
than the EL saturation count rate from NV*^ de¬ 
fects in diamond^®. At higher currents, some 
degradation in the device characteristics is ob¬ 
served which affects the count rate. Therefore, 
no true saturation behavior has yet been ob¬ 
served. 

In contrast, the intensity-power dependence 
in PL shown in Fig. 4b displays a saturation 
type behavior over a count rate similar to that 
observed in EL but with a dependence given by, 
C{P) = CsatP/{P + Psat)- Here, P is the exci¬ 
tation power, C the detector count rate, Psat 
the saturation power and Cgat the saturation 
count rate. Fig. 4b shows a typical intensity- 
power curve with a saturation count rate of 
Csat = 0.9 ± 0.2 MHz and saturation power of 
P^^^ = 880 ± 250 ^W. 

Although the defect emission reaches similar 
count rates in both excitation modes, only the 
PL mode shows the beginnings of saturation be¬ 
haviour. The significantly lower branching ratio 
in EL, indicated by the weak bunching of the 
5 ^^^(r), can therefore yield an explanation for 
the significantly higher count rates without sat¬ 
uration. 

Electrically pulsed single photon emitting 


diode. The single photon emitting diodes allow 
the application of narrow electrical pulses for the 
production of near-deterministic single photons. 
A function generator was used to apply 20 ns 
voltage pulses at 10 MHz in forward bias with a 
peak voltage of 10 V to inject the charge carriers. 
A typical pulsed is shown in Fig. 5a. The 

pulsed emission is indicated by the periodic cor¬ 
relation peaks. The height of the zeroth-order 
peak drops to = 0) ss 0.2, again unam¬ 

biguously indicating single photon emission and 



FIG. 5. Electrically pulsed excitation and 
polarisation dependence, a. Second order auto¬ 
correlation without background correction showing 
g(2)(-r = 0) < 0.2 under 10 MHz pulsed excitation 
where the pulse width and height were 20 ns and 
10 V, respectively, b. Polarisation dependence of an 
electrically driven single emitter. The dashed circle 
represents the background level. 
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in fair agreement with Fig. 3. The non-zero value 
at r = 0 can be accounted for by the background 
EL from the Di center arising from the implan¬ 
tation damage as discussed above. It should be 
noted that this can presumably be further re¬ 
duced with the optimisation of the device geom¬ 
etry. The count rate under these conditions was 
12.2 kHz giving ijdet.Vcoi.Vdef. = 1-2 x 10“^ where 
rjx are the detection, collection and defect pho¬ 
ton generation efficiencies in our system, respec¬ 
tively. An estimation of ijdef. is determined to 
be > 8% with r]det. ~ 53% and rjcoi. ~ 3%. Since 
the device is not yet in the saturation regime 
(Fig. 4), this efficiency is more related to an 
effective recombination efficiency of the defect 
rather than its intrinsic quantum efficiency. 

Excitation repetition rates greater than 
10 MHz decrease the count rate but this may be 
improved once integrated with a cavity or solid 
immersion lens^ since rjcoi. dominates the total 
efficiency of the system. To our knowledge, this 
is the first demonstration of a pulsed room tem¬ 
perature single photon emitting diode. 
Polarization of single centers. A A-half 
wave plate and a polarizer were used to rotate 
and filter the single photon emission. Fig. 5b 
shows the polarisation behavior for a single 
defect. As indicated by the dashed line, the 
emission intensity drops to the background 
level. A study of the polarisation axis of more 
than 30 defects revealed the same fully polarized 
emission along three main polarisation axes 
separated by a 60° angle. This implies that the 
defect emission dipoles are oriented along the 
hexagonal basal plane of the SiC matrix. The 
same behavior was observed in PL in absorption 
and emission. The close similarity of the PL and 
EL polarisation behavior further indicates that 
all defects studied here have the same symmetry 
and therefore suggests the same origin. 

Discussion 

The CMOS compatibility of SiC lends un¬ 
precedented flexibility to further development 
of this electrically driven single photon source. 
This includes further optimisation of the device, 
its integration with optical components such 
as solid immersion lenses, cavities and on-chip 
waveguides and finally, device scale-up. Further 


to this, some flexibility in the integration of other 
defect types, such as the Vsi or the CsiVc, is also 
possible. For instance, defects with a lower ther¬ 
mal stability can be introduced at a later stage 
in the device processing. Optical filtering would 
then be needed to target specific defects. 

Further optimisation of the device design is 
expected to allow higher current densities to be 
generated at the defect site. This may result in 
single photon count rates in the MHz regime. 
The count rate saturation regime (Fig. 3a) 
may then also be accessed to investigate the 
emission dynamics and quantum efficiencies in 
more detail. With these developments and the 
favourable intrinsic properties of these emitters 
such as their full polarisation (Fig. 4b), photo¬ 
stability and room temperature operation, op¬ 
portunities for the application of SiC to quantum 
information processing, quantum cryptography 
and single photon metrology are foreseeable. 

As mentioned above, the spectral charac¬ 
teristics of these sources have not previously 
been reported. We now discuss their possible 
atomic origin. Their characteristic properties 
are their high thermal stability, preferential for¬ 
mation close to the surface, wide variability in 
ZPL and strong polarisation perpendicular to c- 
axis. These common characteristics suggest a 
common origin. Firstly, it is known that stacking 
faults (SFs) often form close to the surface dur¬ 
ing post-growth surface processings^. The vari¬ 
ation in size and type of these SFs could explain 
the large ZPL variabilityS^. However, the polar¬ 
ization of the emitted photons from SFs would 
be emitted parallel to the c-axis because of the 
spatial separation of the electron and hole along 
the c-axis in the excited state^^. This is in con¬ 
trast to the observations made here. Alterna¬ 
tively, the polarisation and high thermal stabil¬ 
ity may be explained by a stable point defect 
with a high Cs^ symmetry. However, a single 
defect type residing on inequivalent sites can¬ 
not account for the observed 200 nm variance 
in ZPL position. We therefore conclude that 
the most likely origin is a combination of these 
two structures, namely, a thermally stable 
symmetry point defect which is near or embed¬ 
ded into a stacking fault where the luminescence 
of the defect involves states split from the con- 
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duction band edge as shown in Fig. 6a. Such a 
model accounts for the spatial location of these 
sources, their polarisation and the fact that their 
density cannot be enhanced by electron irradia¬ 
tion since the SF concentration is intrinsically 
low. Very few defects show a thermal stability 
up to 1600°C in 4H-SiC. One such defect is the 
Di center (Sic shown in Fig. 6b). In our devices, 
luminescence from this center can be clearly ob¬ 
served (for example in Fig. Id) and therefore cer¬ 
tainly has high thermal stability. The Di center 
also has a high Ost, symmetry with polarization 
perpendicular to the c-axis^®. The luminescence 
mechanism can be described as recombination of 
the electron at the split conduction band state 
and a localized hole on the defect site. 

Ab-initio calculations (see Methods) of a Sic 
defect in proximity to a single cubic SF show 
that the ZPL energy is indeed sensitive to the 
defect-SF distance while also maintaining the 
same general recombination behavior found for 
the defect in bulk 4H-SiC. Fig. 6c shows that 
when the defect is placed close to the SF, the 
ZPL energy follows the variation in the conduc¬ 
tion band edge induced by the SF. As the defect- 
SF distance increases the value tends to the bulk 
value. We considered a 2.5-nm wide 3C poly¬ 
type inclusion. This structure yields a triangu¬ 
lar potential well for the conduction band elec¬ 
trons. We find that defects close to but above the 
3C inclusion along the c-axis emit light around 
540 nm. Defects located inside the 3C inclusion 
can emit a continuous range of longer wavelength 
photons from the top side down to the bottom 
side of the inclusion. Defects beneath the 3C 
inclusion show emission around 750 nm. A vari¬ 
ation in defect position near or inside this 3C 
inclusion could then explain the broad variance 
in ZPL position and give rise to the histogram 
displayed in Fig. 2d. 

These calculations also allowed the electronic 
structure and recombination kinetics at the Sic 
defect to be considered in detail (Fig. 6d). The 
highest fully occupied defect level has an e char¬ 
acter in the ground state which lies at about 
0.3 eV above the valence band edge in 4H-SiC, 
and no empty defect state appears in the gap. 
In the excited state the electron occupies an ai 
state split from the conduction band edge and 



FIG. 6. Defect origin, a, Schematic of the band 
edges of a 3C polytype inclusion in 4H-SiC^^. The 
x-axis is parallel to the c-axis with increasing depth 
to the left. The location of a defect state introduced 
by the Sic defect shown in b is symbolized by a red 
dot where the hole is strongly localized upon excita¬ 
tion. c The calculated ZPL wavelength as a function 
of the Sic location around or inside a 2.5 nm wide 3C 
polytype inclusion along, d A schematic of the en¬ 
ergy level diagram of the Sic defect where the ground 
state is ^Ai. Fast (slow) processes are labeled by 
straight (dashed) arrows. PL and EL boxes indicate 
the states participating in the respective processes. 
The state consists of a valence band derived hole 
captured by the defect whereas the ^E and ^E ex¬ 
cited states consist of a hole localized on the defect. 
In these excited states the electron lies on a split con¬ 
duction band state thus the energies of the ^E', ^E 
and ^E levels depend on the position of the conduc¬ 
tion band edge. 


leaves a hole in the strongly localized e state. 
This can form an optically active excited 
state and a lower-energy shelving state. The 
strong luminescence between ^E and ^Ai only 
produces photons with a polarisation perpendic¬ 
ular to the c-axis. Interestingly, the presence of 
this ^E state could explain the photo-dynamics 
we observe in the PL. We also found that the 
e state in the excited state shifts down toward 
the valence band edge compared to the case of 
the ground state. In the EL process the hole 
may first be captured by the valence band de¬ 
rived e state forming a higher-energy ^E' state 
that can decay toward the optically active ^E 






















excited state by spin-orbit interaction mediated 
by phonons. This may explain the four-level- 
like EL process described in Fig. 3. Our working 
model can basically account for all properties of 
the PL and EL single photon sources observed 
in this work. Nevertheless, we do not discount 
other possible candidate extrinsic or intrinsic de¬ 
fects that have a similar electronic structure. 

In conclusion, electrically driven single 
photon sources operating at room temperature 
have been integrated into a range of SiC de¬ 
vices. This is achieved using CMOS compatible 
fabrication techniques, namely photolithogra¬ 
phy, dopant implantation and annealing. The 
high temperature processing resulted in highly 
efficient single photon emitters with a stable 
and high emission rate that exceed that of other 
room temperature electrically driven single 
photon sources. These defect-based sources 
have a high thermal stability, wide variability 
in ZPL wavelength and a strong polarisation. 
The recombination kinetics EL mode showed 
characteristics of a double transition before 
radiative decay while PL model could be char¬ 
acterised by a standard three level energy level 
scheme. A model is proposed which is able to 
account for all of these traits. It consists of 
a Sic in close proximity to a SF-type defect. 
Ab-initio calculations are performed to support 
the validity of this model. The realization of 
bright stable electrically driven single photon 
emission in SiC, in addition to the demonstra¬ 
tion of pulsed excitation of the fully polarized 
single photons, is an important step towards 
on-demand single photons for applications in 
quantum information processing, quantum 
cryptography and single photon metrology. 

Methods 

Device fabrication. The devices were fabricated 
on epilayers grown on 4H-SiC substrates by chem¬ 
ical vapor deposition purchased from CREE. The 
thickness and doping concentration was 4.9 fxm and 
4 X 10^^ cm“^, respectively. Multiple energy A1 im¬ 
plantation was performed into the epilayers while the 
sample was held at 800°C to form the p"*" regions. 
The 6H samples are discussed in the Supplementary 
Information. The implantation energies were chosen 
so as to form a nearly uniform dopant concentration 
depth profile over a 150 nm range. Following implan¬ 


tation an anneal of 1800°C for 10 minutes in an Ar 
ambient was used to activate the implanted dopants, 
repair the implantation damage and also form the 
single photon emitters. The carbon film was removed 
using a dry oxidation process at 800° C following the 
method described in 18. 

A 1100°C pyrogenic oxidation (H2:02 = 1:1) was 
then performed to grow a field oxide. This had the 
effect of bring the SiC surface closer to the peak in 
the implanted dopant profile to allow the formation 
of high quality Ohmic contacts. This oxidation step 
was also found to reduce the PL background. The 
contacts were formed by A1 deposition and subse¬ 
quent sintering at 850°C for 5 minutes. A back con¬ 
tact was formed by A1 deposition without sintering. 
Full electrical characterisation was then performed 
using current-voltage and capacitance-voltage mea¬ 
surements. The PL experiments were performed 
with HPSI 4H-SiC from CREE. Selected samples 
were irradiated with 2 MeV electrons to a fluence 
of 1 X 10^^ cm“^ while the sample temperature was 
kept below 80° C in a N 2 ambient to create a homoge¬ 
neous defect density. Both as-grown and irradiated 
samples were annealed at 1650° C for 10 minutes in 
Ar. Carbonisation of the surface and a subsequent 
1100°C oxidation of these samples were also investi¬ 
gated. 

Confocal fluorescence microscopy. A custom 
confocal system consisting of a HBT interferome¬ 
ter and single-photon-sensitive avalanche photodi¬ 
odes was used to measure the time correlation of 
the fabricated single photon emitters. Photon count¬ 
ing and correlation were carried out using a time- 
correlated single-photon-counting module (PicoHarp 
300, PicoQuant GmbH). The devices were driven 
in forward bias and the EL collected with a xlOO 
infinity-corrected air objective lens with a NA of 0.95. 
The luminescence was collected through a 50 ^m 
multi-mode fibre pinhole. Optical excitation was 
achieved with a 532 nm continuous-wave laser di¬ 
rected through the objective via a dichroic mirror. 
Spectra were collected with a spectrometer (Prince¬ 
ton Instruments Acton 2300i, Acton) with a Peltier- 
cooled CCD. Samples were mounted on a piezoelec¬ 
tric XYZ stage allowing 100 x 100 /rm^ scans. The 
laser line or, during EL experiments, the Di line was 
spectrally removed with a 560-nm long-pass filter. 
The single-photon measurements were performed at 
room temperature. For low temperature measure¬ 
ments the sample was mounted in nitrogen cooled 
cryostat (Linkam THMS600). Instead of a XYZ 
stage, a piezoelectric scanning mirror was used. The 
emission was excited and collected through a cover 
slip corrected 40x air objective with a NA of 0.6. 
Optical excitation was driven by either a 532 nm or 
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a 633 nm continuous-wave laser. The 633 nm laser 
was required to excite emitters with ZPL positions 
at higher wavelengths. 

Modelling. The atomistic simulations were carried 
out using density functional plane wave supercell 
calculations. Our working model is the isolated 
silicon antisite defect (Sic) in bulk and near or 
inside regularly ordered stacking faults, i.e., poly¬ 
type inclusions. We studied the Sic defect at the 
hexagonal site of 4H-SiC in detail. We do not find 
qualitative differences in the ground state properties 
of Sic defect at the cubic or hexagonal sites. The 
Sic defect has a symmetry. In order to eliminate 
the finite size effect of the simulated cell we applied 
a very large supercell consisting of 1536 atoms of 
bulk 4H-SiC in which to embed the defect. We 
applied T-point sampling of the Brillouin-zone that 
is sufficient to provide convergent charge density 
and enabled us to study the degeneracy of the defect 
levels. We applied the Perdew-Burke-Ernzerhof 
(PBE) semilocal DET^® in most of the calculations. 
The PBE functional introduces a band gap error of 
about 0.9 eV in SiC. However, we checked the ground 
state properties of the Sic defect in a 576-atom 
supercell with HSE06 hybrid DET^® which is able 
to reproduce the experimental band gap of 4H-SiC 
(« 3.2 eV) and defect levels in the fundamental 
band gap®°’®^. By comparing the PBE and HSE06 
results we found that a rigid scissor correction on 
the conduction band edge in the PBE calculation 
provided a very good description of the single 
particle levels of the defect state and conduction 
band edge with respect to the calculated valence 
band edge. The value of this scissor correction for 
the zero-phonon line (ZPL) peak of Sic defect is 
0.81 eV. 

The zero-phonon line energies were calculated 
with the constraint DFT (CDET) method where 
the excited state is created by promoting a single 
electron from the occupied state to the empty 
state®^. We also calculated ab-initio the ZPL peak 
of the Sic defect near one cubic stacking fault in 
4H-SiC that may be considered as a half unit of a 
6H polytype inclusion. As a consequence, the band 
gap of the system is reduced by 0.23 eV due to the 
new conduction bands introduced by the stacking 
fault acting as a quantum well for electrons. We 
used a 3016-atom supercell for this study which 
involves five 4H-SiC unit cells and two extra cubic 
bilayers (stacking faults) along the c-direction of the 
lattice. Two stacking faults are needed because of 
the periodic boundary condition along the c-axis. 
These two stacking faults were placed as far away 
from each other as possible (^2.5 nm). We placed 
the Sic defect at a hexagonal site nearest to the 


stacking fault and then three and five bilayers apart. 
We found that the calculated ZPL peaks of these 
defects approach the ZPL peak of Sic defect in 
perfect 4H-SiC as the Sic defects lie farther from 
the stacking fault. The same phenomena is expected 
to arise for larger polytype inclusions such as 3C 
inclusions that should contain a minimum of three 
consecutive cubic bilayers. However, these structures 
are not computationally feasible at the ab-initio 
level. 

For realistically large 3C polytype inclusions in 4H- 
SiC we applied a simplified quasi one-dimensional 
model which adequately describes the electronic 
structure of the defect either near or inside the 
3C polytype inclusions (see 33 for technical and 
other details). The code used to calculate the ZPL 
peaks of polytype inclusions®® was modified so that 
the hole is localized at the defect site whereas the 
electron occupies the lowest-energy empty state in 
the excited state of the defect. The localization 
of the hole is achieved by adding an attractive 
potential for the hole. This attractive potential 
was set to reproduce the experimental ZPL peak of 
the Sic defect in bulk 4H-SiC and then the same 
attractive potential was applied in the polytype 
inclusion model. This methodology was justihed 
from ab-initio calculations on the Sic defect in bulk 
4H-SiC and near to a stacking fault and showed that 
the hole is a deep defect state localized around the 
defect where its level lies at about 0.3 eV above the 
valence band edge. 
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FIG. 1. Atomic structure of hexagonal SiC. Two unit cells of a, 4H-SiC and b, 6H-SiC projected on the 
(1120) plane with lengths of 10.048 A and 15.1248 A, respectively. A projection along <0001> is also shown 
under the 4H-SiC to illustrate the hexagonal configuration of atoms along the basal plane. Kinks along 
the stacking direction, [0001] arising from a vector translation of the bilayers (shaded regions) periodically 
appear and differentiate the polytypes. Point defects of the same type may have different spectral properties 
depending on whether they occupy an inequivalent hexagonal or cubic lattice site denoted by the h and k 
of which there are two in 4H and three in 6H SiC. 


I. MATERIAL CHARACTERISTICS 


Electrically driven single photon sources operating at room temperature have also been inte¬ 
grated into both 4H and 6H SiC devices. These two polytypes are of great technological importance 
and can be grown as virtually dislocation-free 3 inch wafers. The basic structure is a covalently 
sp^ bonded network of alternating Si and C atoms. The various polytypes arise from the different 
ways the SiC bi-layers are stacked along the c-axis as shown in Fig. 1. 

Several batches of devices composed of both 4H and 6H SiC were fabricated using the procedure 
outlined in the methods section of the main text. To predict the profile of implanted species used 
to form the p'’'n or n'’'p junctions we used a Monte Carlo simulation (SRIM code)'^ , the A1 and P 
results are shown in Fig. 2a and 2b, respectively. Multiple energies were used in both cases to form 
a broad constant concentration profile with a peak concentration of 5 x 10^® cm“^. A1 implants of 
110, 75 and 50 keV were performed with fluences of 4.2 x 10^^, 1.3 x 10^*^ and 1.2 x 10^'^ Al/cm^. 
The energies and fluences of the P implants were 140, 90 and 60 keV and 4.2 x 10^^, 1.3 x 10^^ and 
1.2 X 10^^ P/cm^, respectively. During implantation the samples were held at 800°C to minimise 
defect formation and avoid amorphisation. 

Figure 2c and 2d show the forward current characteristics of four different devices with a 1 mm 
diameter top contact and the carrier concentration determined in reverse bias from capacitance 
measurements, respectively. As expected for such SiC diodes, the reverse current of all devices 
tested remained in the sub-nA range for voltages with magnitudes up to 200 V. 
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FIG. 2. Device fabrication and characteristics Monte Carlo simulations of a, the A1 implantation 
profile and b, the P implantation profile after multi-energy implantations, c, The typical forward bias 
current-voltage characteristics for a SiC diodes containing single photon sources. Samples A-D are 6H p^n, 
6H p’*'n, 4H p“''n, and 6H n+p, respectively, d, The carrier concentration versus depletion width in the 
epi-layer of sample C determined with capacitance-voltage measurements to be 9.8 x 10^^ N/cm^. 


II. ELECTRICALLY DRIVEN SINGLE PHOTON SOURCES IN 6H SIC 



FIG. 3. Electrically driven single defects in 6H-SiC. a, Confocal map of a single defect in a 6H SiC 
device. The single defect resides between the floating guard rings which display strong Di emission in this 
device, b, Autocorrelation measurement of the single defect in a. 


As mentioned in the main text, single photon emitters were also integrated into 6H-SiC devices. 
The properties of these emitters were very similar to those found in the 4H-SiC. A typical EL 
confocal map of a single defect in 6H-SiC is shown in Fig. 3a. This particular device is a P 
implanted n+p junction diode. We observed that the floating guard rings were bright with Di 
emission and this is evident as a bright band next to the single defect. The Di center’s PSB 
exhibits a 30 nm red shift with respect to the 4H-SiC devices and so, with a similar hlter to 
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that used for the 4H-SiC measurements, can result in a poorer signal to noise ratio. The high 
recombination rate at the Di center in this device may also explain the general lower emission 
rates of the single photon emitters. 

The measurement of this center is shown in Fig. 3b and proves that it is a single defect. 

The dynamics of defects in 6H-SiC are similar to 4H-SiC, which are presented in the main text. 
This particular defect does not exhibit strong photon bunching, however, the count rates are still 
low compared to the 4H emitters. 

Spectral analysis revealed several different defect types in the spectral region from 600 to 850 nm. 
The spectra were dominated by the very bright Di center emission. Single defect spectra for three 
different defect types are shown in Figure 4d. The spectra are shifted with respect to the 4H 
defects towards longer wavelengths which can be attributed to the same mechanism that causes 
the Di centre emission to shift. As with the 4H devices, PL measurements were not possible due 
to the high doping concentration. 
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FIG. 4. Spectral characteristics of single photon emitters in 4H-SiC. a, Several different but 
reproducible PL spectra excited with a, 532 nm and a, 633 nm. EL spectra for single defects in c, 4H-SiC 
and d, 6H-SiC. 


4 

























III. EL AND PL SPECTROSCOPY OF SINGLE DEFECTS 


A wide range of spectra associated with both electrically and optically driven single photon 
sources is found. A selection of single defect spectra is shown in Fig. 4. A survey of 100 defects 
revealed that many spectral features are located in the 600-700 nm range and are dominated by 
a strong PSB. Fig. 4a and Fig. 4b show PL spectra obtained from single defects in 4H-SiC in 
annealed HPSI SiC under excitation with 532 nm and 633 nm light. The Raman lines are labelled 
TO and LO indicating Raman scattering with the transversal and longitudinal optical phonons. 
The emission spans from the filter edge up to more than 700 nm. The Di line is not exited by the 
laser. Additional EL spectra for 4H-SiC and three selected spectra for single defects in 6H-SiC are 
presented in Fig. 4c and Fig. 4d, respectively. 


IV. POLARIZATION IN PL 




FIG. 5. Polarisation behavior, a, Normalized fluorescence intensity under variation of absorption (black) 
and emission (blue) polarisation filtering for a single defect. A squared sinusoidal curve is fitted and the 
dipole orientation is extracted, b, Left: confocal images under three different, relative laser-polarization 
orientations (red: 30 green: 90 blue: 150 °). Right: merged image. 


To determine the PL dipole orientation in absorption and emission, PL polarization measure¬ 
ments were performed. Figure 5a shows the fluorescence intensity of a single emitter under varying 
excitation polarization (black circles) and its emission polarization under fixed excitation orien- 
tation(blue circles) relative to each other. The squared sinusoidal fit shows a perfect alignment 
between the emission and the absorption dipole. A useful measure is the polarization visibility 
V = {Imax — Imin) / {Imax + Imin), where I is the detector count rate. After subtracting the Raman 
related background (black line), a visibility of Vem = 0.91 was observed for the emission polariza¬ 
tion. The absorption visibility is significantly lower with Vats = 0.59. This behavior was verified 
for most of the investigated emitters. The emission visibility is comparable to that measured for 
EL emitters. 

Figure 5b shows a set of three confocal maps taken under excitation with three different polariza¬ 
tion orientations 60 degrees apart (left). By merging the maps it can be shown that the emitters 
are polarized along three different polarization axes that are 60 degrees apart. 
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FIG. 6. Confocal x — y scan of single photon sources in 4H-SiC as a function of depth. The 

confocal maps are collected with 532 nm excitation approximately from the surface to a depth of almost 
5 /rm. All centers appear in a narrow band close to the surface. 


V. DEFECT LOCATION 

The vertical extent of single photon emitters was determined with a series of vertical x — y 
confocal PL maps taken every 10 nm in depth, z, using optical excitation of the HPSI 4H-SiC 
sample after a 1650°C anneal followed by a 1100°C oxidation. Fig. 3 shows a selection of these 
maps. The superior diffraction limited axial resolution of the PL confocal mode over that of EL 
allowed the approximate depth of the emitters to be determined with a high accuracy. In addition, 
the apparent distribution of defects observed with EL is a convolution of the actual defect density 
with the current paths within the device. In Fig. 6 it can be seen that the defects exist in a narrow 
depth band as their intensity increases and decreases over the range of the map series. When the 
peak intensity of several single emitters versus the depth is fit with a Gaussian distribution after 
corrected for the sample tilt, the emitters have a mean distance from the tilt plane of just ~40 nm. 
This position was common throughout all samples studied and is expected to also be true of the 
pn junction diodes which showed a similar density of alike defects. Therefore, the emitters exist in 
a narrow band close to the surface. 
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